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The near-field scanning optical microscopy images of excitonic wavefunctions in self- assembled 
InAs/GaAs quantum dots are calculated using an empirical pseudopotential method, followed by 
the configuration interaction (CI) treatment of many-particle effects. We show the wavefunctions of 
neutral exciton X° of different polarizations, and compare them to those of the biexciton XX and 
the charged excitons X + and X~ . We further show that the exciton X(Ph — > S e ) transition which 
is forbidden in the far- field photoluminescence has comparable intensities to that of X(Sh — > S e ) 
transition in the near-field photoluminescence . 
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INTRODUCTION 



The 3D confinement and nearly perfect isolation from 
its environment given in self-assemble semiconductor 
quantum dots (QDs) leads to "atom like" electronic 
structure features manifested, among others, by /ieV 
photoluminescence (PL) line widths and a long coher- 
ence time. Such narrow and isolated levels open the 
way for potential applications using as single-photon 
emitters^ and quantum-entangled sources. 2 At the same 
time, this allows fundamental studies of /ieV-scale fine- 
structure splitting^ of exciton lines due to the electron- 
hole exchange interaction, spectral shift due to charged 
excitons£&£ as well as multi-exciton formation and decay. 
However, the spatial resolution of the traditional far-field 
spectroscopy is limited to A/2, where A is the wavelength 
of the probing light, 8 much larger than the typical sizes of 
quantum dots (a few tens to about a hundred nanome- 
ters). Recently, it became possible to exam the inner 
structures of an exciton in the QDs, and map out its spa- 
tial distribution, using the near-field scanning optical mi- 
croscopy (NSOM or SNOM) A±£ In the far field measure- 
ment, the transition intensity (within single-particle ap- 
proximation) is proportional to |p| 2 = |('0 c |r| , 1 ,)| 2 , where 
%j) c and ip v are the conduction and valence band wavefunc- 
tions, respectively. While in the far field measurements, 
the extent of the incident electromagnetic wave is such 
that its magnitude over the entire dot region can be safely 
taken as constant, in the near-field, the electromagnetic 
profile of the tip £(r — ro), has an extent comparable to 
the size of the nanostructure and must be explicitly taken 
into account: 



E(v) = C(r - r )E , 



(1) 



In this case, the transition intensity is modulated by the 
shape of the external field, i.e., proportional to |p(i*o)| 2 = 
|(^c|rf (r - r o)|^)| 2 - The resolution of NSOM can thus 
be much higher than 10 A/2, providing details on how the 
carriers are distributed in the quantum dots. 

In a QD, the electron and hole wavefunctions can 
be approximated by products of envelope functions and 
Bloch functions. In the far field measurement, the inter- 
band transitions obey two kinds of selection rules: one 
related to the Bloch part of the wavefunctions, and one to 
the envelope component of the wavefunctions. For exam- 
ple, from the four possible Sh — > S e excitonic transitions, 
two are bright and two are dark. The dark transitions 
are forbidden because the total angular momentum from 
the electron and hole Bloch functions parts sum up to 
two, and can hence not be carried out by a single photon 
. In contrast, all of the eight Ph — » S e transitions are 
forbidden because a photon does not couple effectively 
to the envelope orbital angular momentum of the exci- 
ton. However, in NSOM, where the electromagnetic field 
is applied inhomogeneously only to a part of the dot, the 
global spatial symmetry is broken, thus altering the se- 
lection rules for the envelope functions and allowing tran- 
sitions that are forbidden in far field measurement^ We 
show here that the Ph — > S e transitions (which are far- 
field forbidden) have comparable intensities in NSOM to 
the Sh — > S e transitions. Furthermore, NSOM provides 
additional far-field allowed information about the Bloch 
part of the wavefunctions, while the far field spectroscopy 
is lacking. For example, by comparing the transition in- 
tensities of different polarizations, one could determine 
the relative phase between \X) and \Y) components of 
the hole wavefunctions at each point. In contrast, in 
the far field measurements, only the total polarization is 
measured and the wavefunction phase information of the 
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Bloch part is lost. 

Recently, Matuda et al^ have measured the near- 
field PL image of naturally occurring GaAs QDs formed 
spontaneously in narrow quantum wells. The enhanced 
NSOM spatial resolution is as high as 30 nm, much 
smaller than the measured QD sizes of ~ 200 nm. Un- 
fortunately, this spatial resolution is still too low for self- 
assembled InAs/GaAs quantum dots. In this paper, we 
predict via atomistic pseudopotential method, the near- 
field PL images of the neutral exciton X°, the biexci- 
ton XX and the charged excitons X + , X~ in the much 
smaller self-assembled InAs/GaAs quantum dots, in an- 
ticipating of future NSOM measurements. 



The a-th neutral exciton wavefunction can be written as, 

N v N c 
v=l c=l 

where N v and N c are the number of valence and con- 
duction states included in the expansion of the exciton 
wavefunctions respectively. The coefficients Cv^} are ob- 
tained by diagonizing the many-particle Hamiltonian H 
in the bassis set {$> vc }i£ 

The far-field optical absorption of the a transition can 
be written as, 



(6) 



II. METHODS 



A. Far- field vs NSOM transition elements 



To obtain the NSOM images, one must first calculate 
calculate the single-particle wavefunctions for electrons 
and holes by solving the Schrodinger equations, 
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"o V 2 + Veff(r) 



(2) 



where V e ff is the effective potential for the electrons in 
the quantum dots. For example, in an effective mass 
approximation (EM A), the effective potential is usually 
taken as two-dimensional harmonic potential, whereas 
in atomistic approaches^ Veff is a superposition of 
screened atomic pseudopotentials. Once one has the elec- 
tron and hole single-particle wavefunctions, it is possi- 
ble to construct the many-particle wavefunctions of the 
excitonic states, e.g., via the configuration interaction 
(CI) method, — by expanding them as the linear combina- 
tion of Slater determinants constructed from the single- 
particle electron and hole wavefunctions. The ground 
state can be described by a single determinant, 

$o(ri,ai, • • • , tat, (Tat) (3) 
= *4[^i(ri,cri) • -i/> v (r v ,<T v ) • • -il) N (r N ,cr N )} , 

where N is the total number of electrons in the system, a 
is the spin index, and A is the anti-symmetrizing opera- 
tor. The Slater determinants of neutral excitons & v c are 
obtained by promoting an electron at the valence state 
ip v to the conduction state ip c from the ground state <I>o, 
i.e., 



where, M^ are the dipole matrix elements between ini- 
tial many-particle state |\E^) and the final state |^/), 

M (fl) = (*/|r|tfi). (7) 

and are the frequency of the external E-field and 
the exciton transition energy, respectively. For a neutral 
exciton, the dipole matrix elements M can be written in 
terms of single-particle orbitals as, 

N v N c 

M (a) = <$o |r|tf x > = EE C t } <^|i# c > , (8) 

v=l c—1 

where, p = (i/j v \t\iIj c ) are the single-particle dipole matrix 
elements between the conduction state iJj c and valence 
state xp v . Similarly, one can calculate M for charged 
excitons X + , X~ and biexciton XX. 

For NSOM, the external E-field is applied locally at 
ro modulated by the electromagnetic profile of the tip 
£(r — ro), which can be taken as a Gaussian function, 
e.g., 



f(r-ro) = c e -l r - r °l 2 / a2 , 



(9) 



where a is the width of the excitation field, much smaller 
than the dimension of the QD. For such local optical 
transition, we have to replace M^ a ^ in Eq. (j6]) by 



N v N c 



M^(r ) =EE^' ftMrtfr - r )|Vc> , (10) 

v=l c=l 

equivalent to replacing in Eq. §8§ the single-particle 
dipole moment p = (^ v \r\i/j c ) by the local single-particle 
dipole moment p(r ) = (^|r£(r - r )|^ c ). 



B. Evaluation of the NSOM transition elements: 
an atomistic approach 



$v,c(ri,<7i, • • • ,r N ,a N ) 

= A[^\{v u (7i ) • ■■i(r CJ (7 c ) • 



(4) In the present work, we calculate the single-particle 

' '0iv( r iV5 & n)] • energy levels and wavefunctions in a atomistic screened 
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FIG. 1: The exciton energies of the lowest neutral exciton X°, 
the charged excitons X + , X~ and the biexciton XX of self- 
assembled InAs/GaAs QD with base 6=27.5 nm and height 
/i=3.5 nm. The right hand side shows the schematic energy 
levels of fine structure splittings of each state. The solid lines 
(Bl, B2, etc.) denote bright (B) states, whereas the dashed 
lines (.Dl, D2, etc.) denote dark (D) states. 



pseudopotential scheme, where the total potential 
V e ff(r) is a superposition of local, screened atomic pseu- 
dopotentials v a (r), plus a nonlocal spin-orbit potential 
V so i.e., V e ff(r) = J2 n ,a M r R n,a) + V so . The atomic 
positions {R U)Q! } are obtained from minimizing the total 
bond-bending and bond-stretching elastic energy using 
the Valence Force Field (VFF) modelJ^^ The atom- 
istic pseudopotentials v a (a=In, Ga, As) are fitted to 
the physically important quantities of bulk InAs and 
GaAs, including band energies, band-offsets, effective 
masses, deformation potentials and alloy bowing param- 
eters, etc. 12 Equation ([2]) is solved by expanding ipfr) 
as the "Linear Combination of Bulk Bands" (LCBB) of 
band index m and wave vector k of material A (= InAs, 
GaAs), strained uniformly to strain (following Ref.[l6l) 
i.e., 



c (A) U W 



,(rK k 



(11) 



We use the lowest conduction bands for electrons, and 
three highest valence bands for holes. Instead of calcu- 
lating the local dipole matrix elements using Eq. (fTT]h 
we project -0(r) onto (strained) Bloch functions of InAs 



at the T point in the Brillouin zone, 



V(r) 



E 



/m(r)ur,m(r) 



(12) 



To simplify the notation, we drop the material index A 
and strain index *e in Eq. ([T2]) . Here, N p is the number of 
the projection bands, which is usually much larger than 
the number of bands used in the LCBB equation. We 
found that N p = 9 provide accurate projections. Here, 
/m(x) is the envelope function for the ra-th band. The 
envelope functions are slow varying functions, compared 
to the Bloch functions. We can therefore separate the 
Bloch functions and envelope functions by dividing the 
volume into small regions, e.g. the 8-atom unit cells. In 
each unit cell, the envelope functions / m (x) and £(r) can 
be treated as constants, i.e., 



dr^*(r)V> c (r)r£(r - r ) 



E 



E /^ ) (r,)/i c , ) (r i )d m ,, 



m,m' 



(13) 



f (r* - r ) 



where, is the position of i-th cell and 

d m,m' = ^7 j dr ^r 5m (r)r^ r ,m'(r) , (14) 

is the dipole matrix element between bulk bands. In a pe- 
riodic system, it is more convenient to use an alternative 
formula to calculate the dipole, 17 
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imUJ rn r n > 



-(ur,m\p\ur,m>) • 



(15) 



If we assume high spatial NSOM resolution, i.e., £(r^ 
r o) = S(ri — i*o), Eq. (|T4|) can be further simplified as, 



P(ro)« E rt v) (ro)f£(ro)d r , 



(16) 



Combining Eq. ((TDj) and Eq. (p~6|) . we see that the shape 
of NSOM transition is the sum of the (multi-band) en- 
velope functions of the exciton weighted by the (single- 
particle) dipole matrix elements, i.e., 



N v N c 



v=l c—1 



(17) 



Previously, the NSOM images have been calculated via 
single-band effective mass approximation (EMA) i n i 18 i 19 
In this approximation, the dipole moment d is directly 
taken from bulk materials or from quantum wells, 19 while 
in our case, it is calculated for nine bands in Eq. ([T5]) and 
is specific to the nanostructure considered. In the single- 
band EMA, only the envelope part of the information is 
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kept in the calculation of NSOM images, whereas detailed 
information associated with the Bloch functions, such as 
polarizations is lost. In the more advanced k-p methods, 
one includes the heavy- hole, light-hole, and spin-orbit 
bands at T point of the Brillouin zone in the calculations. 
For the direct gap materials, these approximations give 
qualitative right results, but for indirect gap materials, 
such Si quantum dots, one need to include many more 
bands. 



III. RESULTS AND DISCUSSION 



A. The excitonic spectrum of the InAs/GaAs QD 



We performed calculations on a lens-shaped 
InAs/GaAs quantum dot with 27.5 nm base diam- 
eter and 3.5 nm height. The calculated transition 
energies of the exciton X°, the charged excitons X~ 
and X + and the biexciton XX of Sh — > S e transitions 
are given in Fig. [TJ The neutral X° has the largest 
transition energy of 981 meV, while charged excitons 
and biexcitons have 1-3 meV red shifts relative to the 
neutral exciton, due to correlation effects, 20 in agree- 
ment with experiments^ For the neutral exciton, the 
exciton line splits into four lines shown on the righthand 
side of Fig. [TJ as a consequence of the asymmetric 
electron- hole exchange energies The two low lying 
lines labeled Dl, D2 are "dark states" (forbidden due 
to the angular momentum selection rule of the Bloch 
part of the electron/hole states), whereas the upper two 
states (51, B2) are bright states. The calculated energy 
splitting between the bright and dark states is about 210 
fieV. The energy splitting between two dark states Dl, 
D2 is about 1 /ieV, whereas the splitting between two 
bright states Bl, B2 is about 10 /ieV, in agreement with 
previous calculations *2£i The low energy bright state Bl 
is polarized along the [110] direction, whereas the high 
energy state B2 is polarized along the [110] direction. 
Although the biexciton XX states themselves have no 
fine structures, the XX to X transition exhibit fine 
structures. The charged excitons X + and X~ have no 
fine structures. 

The energy of neutral exciton Ph — > S e transition is 
calculated here 995 meV, also shown in Fig. [TJ Like the 
Sh — > S e transition, the Ph — > S e transition also split 
into four fine structures. The two lower energy states 
are "forbidden" by the Bloch function angular momen- 
tum selection rule, whereas the two high energy states 
Bl and B2 are allowed by the Bloch part. In the far 
field measurement, the Ph — > S e transition is forbidden 
for the lens-shaped QDs, due to the envelope function 
angular momentum selection rule (i.e., the overlap in- 
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FIG. 2: (Color online) The NSOM of excitonic S h -> S e tran- 
sitions in self- assembled InAs/GaAs QD with 6=27.5 nm and 
/i=3.5 nm, calculated from Eq. (|17|) . The two dark excitons 
Dl, D2 are plotted in (a) and (b) respectively, whereas the 
two bright excitons Bl, B2 are plotted in (c) and (d). The 
unit of the color-bar is arbitrary but equal for all sub-figures. 
In each sub-figure, there are six panels showing the transi- 
tion intensities of different polarizations. In the first row, the 
polarizations are along the [100], [010] and [001] direction re- 
spectively, and in the second row, the polarizations are along 
the [110] and [110] directions. The last panels show the sum 
of transition intensities of different polarizations. We magnify 
the transition intensities by the factors on the upper-right 
corner of each panel. The white circles (27.5 nm in diameter) 
show the boundaries the QDs. 
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FIG. 3: (Color online) The NSOM of excitonic P h -> S e 
transitions (see Fig. [T} in self- assembled InAs/GaAs QD with 
6=27.5 nm and /i=3.5 nm, calculated from Eq. ([17)1 for (a) 
Bl and (b) B2 states. The unit of the color-bar is arbitrary 
but equal for all sub-figures. In each sub-figure, there are six 
panels showing the transition intensities of different polariza- 
tions. In the first row, the polarizations are along the [100], 
[010] and [001] direction respectively, and in the second row, 
the polarizations are along the [110] and [110] directions. The 
last panels show the sum of transition intensities of different 
polarizations. We magnify the transition intensities by the 
factors on the upper-right corner of each panel. The white 
circles (27.5 nm in diameter) show the boundaries the QDs. 



tegrals between Ph and S e envelope functions are very 
small). In practice, this transition is partially allowed, 22 
but nevertheless very weak. However, unlike the selection 
rules for the Bloch functions, the envelope function an- 
gular momentum selection rule does not apply to NSOM, 
which involves only the local transition and does not feel 
the global spatial symmetry^ i.e., the overlap of Ph and 
S e orbit als are not zero at given point, and Ph — > S e 
transitions have comparable intensity to the Sh — > S e 
transitions. This will be further discussed in Sec lIII CI 



B. NSOM Sh — ► S e transition for neutral exciton in 
the InAs/GaAs QD 



The NSOM images of excitonic Sh — > S e transition of 
X° are shown in Fig. [2j We show the images of all fine 
structures components Bl, B2, Dl, D2 in Figs. E^a), 



[2(b), [2(c) and [51(d) respectively. The color-bars indicate 
the magnitudes of transition intensities. We use arbitrary 
but equal units for the transition intensity in all sub- 
figures. Each sub-figure is composed of six panels show- 
ing the transition intensity of different polarizations. The 
first row in each sub-figure shows the transition of polar- 
ization along the [100], [010] and [001] directions, whereas 
the second row of each sub-figure shows the transition of 
polarization along the [110] and [110] directions. The fi- 
nal panel is the total transition intensity, being the sum 
of intensities over polarization along the [100], [010] and 
[001] directions. (The sum of intensities over the [100], 
[010] polarization equals the sum over the [110] and [110] 
polarization.) We magnified the transition intensity by 
the factor shown on the upper right corner of each panel, 
so the transition intensities in all sub-figures can be in 
similar scale. 

Considering first the total transition intensities, we no- 
tice that the transitions originate from a relatively small 
region around the center of the dot. The transition in- 
tensities of dark states Dl [Fig. [2(a)] and D2 [Fig. [2(b)] 
are about 2-3 orders of magnitude lower than the those 
of bright exciton states Bl [Fig. [2(c)] and B2 [Fig. [2(d)], 
because they are "forbidden" due to angular momentum 
selection rule for the heavy-hole to conduction band tran- 
sitions. Therefore the shapes of Dl and D2 are deter- 
mined by the shape of higher bands [see Eq. ([16])]. In 
contrast, the bright states Bl and B2 have S- like shape 
as expected. However, the NSOM of Bl and B2 are 
not perfectly round, even though the lens-shaped QD we 
choose has cylindrical symmetry: we see that Bl is elon- 
gated along the [110] direction, whereas B2 is elongated 
along the [010] direction. 

We next discuss the NSOM images of different polar- 
izations, focusing on the bright states. For far field exci- 
tation, the low energy bright exciton Bl is calculated to 
be polarized along the [110] direction, whereas the high 
energy bright exciton B2 is polarized along the [110] di- 
rection. This is clearly seen in the NSOM plots of Fig. 
[2j by comparing the [110] polarized plots to the [110] po- 
larized plot. However, NSOM can provide much more 
information than the usual far field measurement: The 
dipole moments have the following relations, P[no] = 

(P[ioo] +P[oio])/v / 2 and p [li0 ] = (P[ioo] - P[oio])/>/2. 
Therefore, by comparing the NSOM images of different 
polarization, we can obtain the relative phases between 
P[ioo] an d P[oio] a t eacn spatial point. For the bright 
states of Sh — >• S e transitions, the relations are simple, 
i.e., P[ioo] an d P[oio] have almost fixed phases ir/2 every- 
where inside the QD. As we know, the P[ioo] correspond- 
ing to the \X) component in the hole Bloch functions, 
whereas P[oio] corresponding to the \Y) component. We 
therefore obtain additional information about the char- 
acter of the hole Bloch functions as a function of position 
inside the QDs. 
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Figure [2] also show the NSOM image of polarization 
along the [001] direction. For the bright exciton, the 
[001] components are significantly weaker than the com- 
ponents polarized in (OOl)-plane, and may not be de- 
tectable in practice. However, they are of theoretical 
interest. The polarization along the [001] direction come 
from the emission of light hole. Therefore, the [001] im- 
age shows how the light hole components distributed in 
the QDs. 



C. NSOM Ph — > S e transition for neutral exciton in 
the InAs/GaAs QD 



Figure [3] depicts the Ph — > S e transitions, showing the 
two bright states (Bl, B2) out of the four fine-structure 
split states. Similar to Fig. O we show the NSOM images 
with polarization along the [100], [010], [001] directions in 
the first row, and those with the [110], [110] polarization 
and the total transition intensities in the second row. 
We see that the Ph — > S e transitions have comparable 
intensities to those of the Sh — > S e transitions. The shape 
of total transition intensity of the Bl state somewhat 
differs from the intensity from the B2 state. However, 
both states have ring like structures, and both states have 
two peaks in the [110] direction. 

Each of the polarization-resolved plots looks like 
orbital-p functions, although the transition peaks are ori- 
ented in different directions. For example, for the Bl 
state, the peaks of the [100] ([010]) polarized light are 
oriented along the [100] ([010]) direction, whereas the 
peaks of [110] ([110]) polarized light are oriented in [110] 
([110]) direction. As discussed in Sec lIII B\ we can ob- 
tain the phase relations of \X) and |Y) components in 
hole wavefunctions. We see that unlike the previous case 
involving an Sh wavefunction, in the Bloch part of the 
Ph wavefunctions, \X) and \Y) components have differ- 
ent phases at each point of the QDs. Similar results are 
found for the B2 state. 



D. Comparison of NSOM of X°, X + , X~ and XX 
excitons in the InAs/GaAs QD 



In Figs. H (b), (c), (d) we depict the NSOM images 
of the charged exciton X + , X~ and the biexciton XX, 
respectively, and compare them to the X° Sh — > S e tran- 
sition from Fig. UK a). We also plot the NSOM transition 
intensities along the [110] direction in Fig. [5] for X°, X + , 
X~ and XX. The transition intensities of X°, X + , X~ 
and XX along [110] direction are very similar to those 
along the [110] direction, and therefore are not plotted. 
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FIG. 4: (Color online) The total NSOM images (sum over 
all fine structures) of excitonic Sh — ► S e transitions (see Fig. 
QJ for (a) X°, (b) X + , (c) X~ , and (d) XX in self- assembled 
InAs/GaAs QD with 5=27.5 nm and /i=3.5 nm, calculated 
from Eq. ([IT]) . The unit of the color-bar is arbitrary but equal 
for all sub-figures. In each sub-figure, there are six panels 
showing the transition intensities of different polarizations. 
In the first row, the polarizations are along the [100], [010] 
and [001] direction respectively, and in the second row, the 
polarizations are along the [110] and [110] directions. The 
last panels show the sum of transition intensities of different 
polarizations. We magnify the transition intensities by the 
factors on the upper-right corner of each panel. The white 
circles (27.5 nm in diameter) show the boundaries the QDs. 
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FIG. 5: (Color online) The NSOM transition intensities of X° 
(black solid line), X + (red dashed), X~ (green dashed) and 
XX (blue dashed) along the [110] direction in self- assembled 
InAs/GaAs QD with 5=27.5 nm and h—3.5 nm, calculated 
from Eq. (|17|) . The origin is at the center of the dot. 



We see that X has the strongest transition intensity, 
followed by the X°, XX and X+. Even though, the 
charged excitons and the biexciton have 1-3 meV red 
shifts relative to X°, they have very similar shapes and 
sizes to X°. The only noticeable difference being that 
X + has a slightly lager size than the other transitions. 
In Ref. 9, it was found experimentally and confirmed the- 
oretically by single-band effective mass calculation s 18 ! 19 
that the radius of the biexciton NSOM signal is smaller 
than the radius of X°. In this calculation, we did not 
find a significant difference in size and shape between X° 
and XX. However, the QDs used in the experiments 
(created by composition fluctuations in a quantum well) 
are much larger (~ 200 nm in diameter) than the self- 
assembled QDs we studied here, therefore the correlation 
effects are stronger. On the other hand, we ignored the 
continuum states in the CI calculations, which may also 
underestimate the correlations effects to the wavefunc- 



IV. SUMMARY 



We calculate the NSOM images of exciton complexes in 
self- assembled InAs/GaAs QD using an atomistic pseu- 
dopotential method followed by a configuration interac- 
tion treatment of the correlations. The NSOM images 
reveal the spatial structure of exciton and give informa- 
tions about the Bloch function character of the wavefunc- 
tion, which remain unobserved in conventional far-field 
measurements. We show that the — » S e transition, 
which is forbidden in the far field PL measurement, has 
comparable transition intensity to that of the Sh — ► S e 
transition in NSOM measurements. We also calculate the 
NSOM image of charged exciton X + , X - , and biexciton 
XX, and compare them with the image of X°. We found 
that the images obtained for charged excitions and the 
biexcition show very similar features to the ones obtained 
from the exciton decay. 
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